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In  a  previous  study  we  have  developed  a  general  theory  of  chemisorption  on 
certain  stepped  surfaces;  namely  those  for  which  a  lattice  gas  in  thermo¬ 
dynamic  equilibrium  is  an  appropriate  model.  We  showed  there  that  a  comparison 
of  LEED  (Low  Energy  Electron  Diffraction)  results  for  flat  and  stepped  substrates 
can  determine  the  change  in  adsorption  energy  at  terrace  edge  sites.  In  the 
present  work  we  demonstrate  this  via  detailed  Monte  Carlo  calculations  for 
0  on  a  certain  stepped  W(110)  surface.  We  use  oxygen  adatom-adatom  (AA)  in¬ 
teraction  energies  previously  determined  by  experimental  and  theoretical  studies 
on  flat  W(110).  By  comparison  with  experimental  results  for  the  stepped  surface 
system  we  find  that  the  0  binding  energy  is  less  strong  at  either  terrace  edge 
for  this  particular  surface  and  defect. 

1 

I.  Introduction 

In  a  previous  study1  (referred  to  as  .  herein)  we  have  presented  a  general 
theory  of  LEED  scattering  from  overlayers  on  stepped  surfaces,  including  the 
effects  of  statistical  disorder  in  the  overlayer.  Systems  that  may  be  modeled 
by  a  lattice  gas  in  thermal  equilibrium  were  considered.  We  showed  that  the 
behavior  of  the  LEED  intensity  as  a  function  of  coverage  0  is  strongly  influenced 
by  the  sign  of  the  change  in  adsorption  energy  at  terrace  edge  sites.  We  also 
considered  LEED  spot  splitting  effects  and  pointed  out  that  qualitatively  new 
LEED  features  may  arise  at  low  (or  high)  0  values,  due  to  differences  in  adatom- 
adatom  (AA)  interaction  energies  at  terrace  edge  sites  from  their  value  else¬ 
where  on  the  terrace. 

♦Supported  in  part  by  the  Office  of  Naval  Research  and  a  Stauffer  Chemical 
Company  Grant  of  Research  Corporation. 

♦♦Present  address:  Department  of  Mathematics,  SUNY-Buffalo,  Buffalo,  N.Y. 
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The  present  work  applies  the  theory  developed  in  I  to  existing  experimental 
o 

results  for  0  adsorbed  on  a  certain  stepped  W(110)  surface.  In  Section  II  we 
briefly  review  the  theory  presented  in  I.  In  Section  III  we  describe  the  ex¬ 
perimental  data  for  0  on  stepped  W(110).  Integrated  LEED  intensity  as  a  function 
of  coverage  0  and  spot  splitting  results  for  the  two  adlayer  domains  are  described. 
Section  IV  defines  our  theoretical  model.  A  Monte  Carlo  program  was  employed  to 
calculate  the  statistical  averages  of  interest.  As  input  we  used  AA  energies  for 
0/W(110)  derived  by  comparison  of  experimental  studies  on  the  flat  surface  and 
theoretical  lattice  gas  models.  Section  V  contains  the  results  of  our  calculations. 
Using  the  signs  of  the  changes  in  adsorption  energy  at  terrace  edge  sites  as 
unknown  parameters,  we  find  that  the  observed  LEED  features  are  reproduced  most 
accurately  when  the  0  adsorption  is  less  strong  at  either  terrace  edge  site 
than  on  the  flat  surface.  Section  VI  contains  some  suggestions  for  future  work. 
Perhaps  most  exciting  is  the  prospect  -  already  suggested  by  our  results  in  I  - 
that,  very  generally,  one  can  determine  the  sign  of  the  change  in  adsorption 
energy  at  terrace  edge  sites  by  a  simple  comparison  of  appropriate  overlayer 
LEED  results  for  the  flat  and  stepped  surfaces  without  detailed  modeling. 

II.  Theory  of  LEED  Scattering  from  Overlayers  on  Stepped  Surfaces 

In  this  section  we  review  the  main  points  of  the  theory  developed  in  I. 

There  we  assumed  the  stepped  surface  was  composed  of  identical  terraces,  with 
equivalent  points  on  neighboring  terraces  separated  by  a  fixed  vector  £,  and 
that  N$  terraces  were  within  the  coherence  length  £  of  the  LEED  apparatus.  We 
assumed  the  adatoms  on  a  given  terrace  sit  at  N  fixed  sites  R^,  with  occupation 
variable  n^  =  0(1)  for  an  unoccupied  (occupied)  site.  If  the  adatoms  are  in 
thermal  equilibrium,  and  the  AA  forces  of  not  too  long  range,  it  is  reasonable 
to  assume  that  the  overlayer  on  each  terrace  is  statistically  independent  of  that  on 
the  others.  Another  important  point  is  that  for  strong  chemisorption  the  change 


"1 


V 


1 

6e  in  adsorption  energy  at  terrace  edge  sites  will  generally  be  such  that 
|6e|  »  kT.  Hence  edge  sites  are  either  empty  (Se  >  0)  or  occupied  (6e  <  0) 

* 

given  a  sufficient  value  of  9. 

i 

The  result  of  all  this  was  that  in  the  kinematic  (single-scattering)  ap¬ 
proximation  the  LEED  intensity  for  scattering  vector  k^  is 

I(ls)  =  N,(-elf<k>lS  -  U}|<y(fcV>r  L\) 

Here  ^ 

s  ,  6) 
s  siv?(Nsx/i)  /  sin*  (jr/i)  . 

N  is  the  number  of  adsorption  sites  on  a  terrace  and  the  angular  brackets  denote 
a  statistical  average  on  a  single  terrace.  The  first  term  in  Eq (1 )  includes 
effects  of  statistical  disorder  on  a  single  terrace,  while  the  second  term 
includes  effects  of  interterrace  interference  via  the  function  h.  Thus  for  a 
given  stepped  chemisorption  system  and  set  of  LEED  scattering  parameters,  the 
intensity  I  is  determined  by  the  temperature  T,  coverage  0,  AA  interactions, 
and  signs  of  the  changes  6e  in  adsorption  energy  at  terrace  edge  sites.  All  of 
these  quantities  except  <5e  may  be  measured  or  determined  by  calculations  for 
the  corresponding  flat  surface  adsorption  system.  Thus  the  sign  of  6e  can  be 
found  by  studying  experimental  results  for  the  stepped  surface.  This  is  shown 
by  means  of  a  model  calculation  at  T  =  0  in  I  and  also  demonstrated  with  a 
(finite  temperature)  Monte  Carlo  calculation  for  0/stepped  W(110)  here. 

In  I  we  also  consider  several  features  of  the  LEED  scattering  in  more 
detail.  In  particular,  we  show  in  a  model  calculation  that  the  overlayer  beam 
intensity  I(k^'),  where  Is  an  appropriately  chosen  scattering  vector,  has 
a  maximum  at  a  coverage  0m  that  is  characteristic  of  the  change  in  adsorption 
energies  6e.  This  result  is  only  due  to  the  change  in  effective  number  of 
adsorption  sites  available  on  the  stepped  vs.  flat  surface.  Thus  one  would 
expect  a  similar  trend  in  the  maximum  of  the  integrated  intensity  -  vs  9,  where 


In  Eq(4),  and  J)  are  unit  mesh  vectors.  •<  is  an  experimentally  measurable 

quantity.  The  reciprocal  space  integration  area  ou is  determined  by  instrumental 

parameters  and  beam  energy.  This  expectation  is  confirmed  by  the  results  in  Section  V. 

2 

In  I  we  go  on  to  consider  the  relative  size  of  the  contribution  of  |<p>| 

2  1 

and  < | p |  >  to  >  (these  terms  enter  via  the  expression  for  I  given  In  Eq(l)). 

We  argue  that  the  former  often  tends  to  be  smaller.  We  also  consider  the  splitting 
of  LEED  spots,  which  is  due  to  interterrace  interference  via  the  function  h  in 
Eq (1 ) .  We  show  by  means  of  general  arguments  and  specific  model  calculations 
that  the  splitting  strength  depends  on  coverage,  the  signs  of  the  fie,  and  whether 
the  terrace  width  is  odd  or  even. 

In  what  follows  we  consider  the  chemisorption  systems  0/W(110)  and  0  on  a 
certain  stepped  W(110)  surface.  We  show  that  applying  the  general  theory  of 
I,  described  above,  to  experimental  results  for  these  systems  leads  to  the  con¬ 
clusion  that  the  0  adsorption  energy  is  less  strong  (6e>0)  for  either  terrace 
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edge  in  this  system.  Similar  information  has  been  obtained  in  field  emission 
studies  on  other  systems.  There  seems  to  be  no  general  rule  as  to  whether  a 
given  step  is  more  or  less  adsorptive  for  a  given  adsorbate. 

In  applying  our  theory  we  make  use  of  the  AA  interaction  energies  derived 

4 

by  Ching  et  al  in  a  lattice  gas  treatment  of  the  flat  surface  chemisorption  system. 
These  parameters  are  employed  as  input  in  a  statistical  mechanical  calculation 
of  LEED  scattering  on  the  stepped  surface  chemisorption  system  with  the  fie  as 
variables.  The  results  of  this  calculation  are  compared  with  experiment.  We 
find  agreement  only  for  fie>0  on  both  terrace  edges. 

III.  Experimental  Results 

c 

The  system  0/W (110)  has  been  extensively  studied  .  For  0  <  0.5,  and 

T<Tr(0)  the  adlayer  shows  (2x1)  ordering.  For  T  >  T  it  is  In  a  disordered 
c  c 


r 


state.  The  phase  boundary  T  (6)  Is  observed  to  be  approximately  constant, 

C 

T  =  450°K  for  .1  <  0  <  0.35  and  then  jump  abruptly  to  T  s  720°K  for 

L  —  w 

0.35  <  0  £  .5.  AA  (adatom-adatom)  Interaction  energies  have  been  derived  by 
several  authors  by  comparing  experimental  and  theoretical  T  results  at  a  few 

v 

0  values.  We  use  the  values  determined  by  Ching  et  al  which  are  illustrated 
in  Fig.  1. 

Note  that  in  this  system  there  are  two  degenerate  (2X1)  structures  which 
are  completely  equivalent  in  the  absence  of  defects.  Hence  the  (1/2, 1/2)  and 
(1/2, 1/2)  LEED  beams  have  the  same  behavior  as  functions  of  0  and  T. 

2 

An  experiment  for  0  on  stepped  W(110)  has  been  performed  by  Engel  et  al. 

The  W  crystals  were  cut  so  that  the  average  terrace  width  D  on  one  was  D=10 
adsorption  sites  and  the  other  D=24.  In  the  nomenclature  of  Lang  et  al.®  these 
can  be  indexed  as  W-(S)  [  10(1 1 0)x(0l 1 )]  and  W-(S)  [24 (llO)x(Oll )].  Thus  the 
terrace  edges  are  in  the  [Tl 0]  direction,  as  illustrated  in  Fig.  2.  This  means 
that  the  two  (2x1)  overlayer  domains  become  crystallographically  inequivalent. 
Domain  I  refers  to  that  with  the  [0,1]  direction  parallel  to  the  terrace  edge, 
domain  II  to  the  other,  as  shown  in  Fig.  2.  The  (1/2, 1/2)  ((1/2, 1/2))  LEED 
beam  corresponding  to  domain  I ( II )  was  monitored  in  this  experiment. 

In  the  calculations  reported  here,  we  considered  the  D=10  case  only  to 
minimize  computer  time.  However  the  experimental  results  at  D=24  are  completely 
consistent  with  our  interpretation.  At  D=10,  experiments  were  done  at  three 
temperatures:  T=300°K,  460°K  and  1000°K.  We  calculated  at  460°K  only. 

There  Is  evidence7  that  0  diffusion  rates  on  this  surface  fall  considerably 
below  300°,  so  that  thermodynamic  equilibrium  may  not  be  achieved.  For 
T  >  900°K  the  steps  on  the  clean  surface  begin  to  move2.  (We  note,  however, 
that  the  results  reported2  for  T  =  460°K  and  900°K  are  in  qualitative  agreement 
with  our  interpretation.) 
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The  experimental  results  we  will  deal  with  Include  (a)  the  integrated 
LEED  intensity  •  and  (b)  the  LEEO  overlayer  spot  shapes.  Under  (a)  it  was 
observed  that  while  the  integrated  intensities  /  j  (for  domain  I)  and  '.jj 
(for  domain  II)  are  roughly  equal  at  low  coverage  (0<O.l), Tj  is  considerably 
larger  than  tjj  for  larger  values  of  0,  as  one  would  expect  from  the  terrace 
orientation.  Further,  the  maximum  value  of  -j  is  reached  not  at  0  =  0.5, 
which  would  correspond  to  a  complete  (2x1)  overlayer  on  a  flat  surface,  but  at 
0  =  0.32.  (for  D=24  the  corresponding  maximum  is  at  0  =  0.44).  has  a 
minimum  at  a  coverage  slightly  below  the  maximum  in  in  both  cases.  These 
results  are  illustrated  in  Fig.  3. 

Under  (b),  the  (1/2, 1/2)  (domain  I)  beam  shows  a  splitting  in  the  [T  1  2] 
direction  at  characteristic  energies,  from  which  terrace  width  and  step  heights 
were  deduced.  Also  a  minimum  in  intensity  at  the  normal  petition  of  the  (1/2, 1/2) 
reflex  was  seen.  The  (T/2,1/2)  beam  (domain  II)  showed  no  splitting.  All 
overlayer  features  were  streaked  in  the  {.T  12]  direction,  which  is  consistent 
with  a  distribution  of  terrace  widths  and  the  clean  surface  LEED  pattern. 


IV.  Theoretical  Model:  Details 

We  have  calculated  LEED  beam  intensities  I(k)  and  integrated  intensities 

* 

using  the  theory  described  in  Section  II  in  a  model  appropriate  to  the 
0/stepped  W(110)  results  discussed  above.  The  Monte  Carlo  technique  was  used. 
All  calculations  were  done  at  T  -  0.04  eV  on  a  single  terrace  of  size  Dx30, 
with  periodic  boundary  conditions  in  the  if  11]  direction.  Our  program  was 

O 

similar  to  one  described  by  Landau  .  The  grand  canonical  ensemble  was 
employed,  and  we  used  6  randomly  chosen  Initial  configurations,  calculated  240 
configurations  in  each  stream  and  kept  140  for  averaging.  The  results  reported 
are  quite  independent  of  the  Initial  configuration  -  full  or  empty  initial 
configurations  changed  none  of  the  quantities  significantly. 
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According  to  the  general  theory  given  above,  there  are  three  possibilities 
for  the  change  in  adsorption  energy  6e  at  the  terrace  edges:  both  <Se>0  (+  +), 
both  6e<0  ( — )  and  one  of  each  (+  -).  We  set  |6e|  =  leV  in  all  cases.  To 
obtain  <p(k)>  we  Fourier  transformed  <n.>  directly.  For  (+  +)  amd  (-  -)  these 
averages  were  symmetrized  after  calculation.  The  integration  region  for  Eq(4) 

p 

was  not  reported  by  Engel  et  al.  For  domain  I  (II)  we  assumed  it  to  be  an  equal - 

sided  parallelogram  centered  on  the  (1/2, 1/2)  ((1/2, 1/2))  beam  with  an  area  of  4%  of 

the  (flat  surface)  Brillouin  zone.  In  doing  this  h  in  Eq(l)  was  set  equal  to  its 

maximum  at  the  adlayer  spot  center.  Since  k*£  was  not  reported,  we  set  it  equal  to 

n  at  the  adlayer  spot  center  and  took  £  '{  a  in  the  integration  in  Eq(4).  These 

assumptions  amount  to  dividing  the  [T  12J  direction  equally  into  regions  centered 

on  the  peaks  of  h.  The  integral  of  h  is  equal  to  the  integral  of  N$  over  such 

a  region,  and  we  found  the  contribution  of  | <p> (  to\rto  be  small.  This  confirms 

2  2 

the  general  argument  advanced  in  I  (see  Section  II).  <|p|  >  and  |<p>|  are 
similar  smooth  functions  of  9  so  a  change  of  integration  region  in  vshould  not 
alter  our  main  conclusions,  which  are  based  only  on  trends  in  this  quantity. 

However,  the  precise  value  of  the  coverage  0  at  the  peak  in  r.j  might  be  altered 
slightly. 

We  note  that  a  very  interesting  renormalization  group  treatment  of  finite 

g 

size  effects  on  overlayer  thermodynamics  has  been  presented  by  Berker  and  Ostlund  . 
This  gives  an  alternate  means  of  calculating  some  of  the  quantities  determined 
with  a  Monte  Carlo  treatment  here. 

V.  Theoretical  Model:  Results 

First  we  discuss  the  integrated  LEED  intensity. 

Results  for  j  vs.  9  are  shown  in  Fig.  4  for  six  cases:  terrace  width 
D=9  or  10  and  terrace  edge  adsorption  energy  changes  (+  +),  ( — ),  and  (+  -). 

We  chose  both  an  odd  and  even  value  of  D  near  the  observed  average  to  check 


1 

for  odd-even  effects,  which  can  be  significant,  as  pointed  out  in  Section  II. 

The  (-  -)  cases  show  no  peak  in  ^  j  for  0  <  0.5.  For  D=9,  (+  -)  there  is  a  peak 
slightly  below  e  =  0.5.  However -  j  for  D=10,  (+  -)  is  increasing  in  this  6 
region  and  an  average  of  the  two  shows  little  or  no  peak.  This  average  is 
appropriate  to  the  real  surface,  as  remarked  above.  Clearly  the  (+  +)  cases 
best  reproduce  the  most  striking  feature  in  j,  a  maximum  for  0  <  0.5.  We 
regard  this  as  the  most  convincing  piece  of  evidence  for  weaker  adsorption  at 
the  terrace  edges.  The  apparent  reason  for  it  is  also  quite  intuitive  -  excluding 

p 

U  of  the  adsorption  sites  means  that  a  complete  (2x1)  adlayer  is  formed  at  a 
lower  value  of  0.  (This  remark  is  elaborated  on  below,  also  see  I).  The 
experimental  peak  comes  at  a  lower  0  value  than  those  in  Fig.  4.  In  the 
context  of  our  model,  there  are  several  possible  reasons  for  this  discrepancy. 

(i)  The  value  of  the  Auger  signal  defining  0  *  0.5  was  taken  to  be  that  at 
which  the  LEED  signal  for  the  (2x1)  domain  on  the  flat  surface  was  a  maximum. 

This  calibration  procedure  is  somewhat  sensitive  to  the  defect  concentration  on 
the  flat  surface,  and  also  to  instrumental  effects.  ( i 1 ) .  The  value  of  the 
LEED  signal  integration  area  a.  on  the  stepped  surface  was  not  reported,  and 

i 

the  peak  in  ^  may  depend  slightly  on  this  quantity,  as  mentioned  above. 

(iii).  The  stepped  surface  used  had  a  distribution  of  step  widths  D,  and  the 

.  2 

shift  in  the  peak  in  j,  depending  mainly  on  ^  (see  below),  is  larger  for 
smaller  D  values,  (iv).  The  calculation  presented  here  is  for  a  single  value 
of  D  only.  (v).  A  final  possibility  is  that  the  adsorption  energy  decreases 

at  terrace  edge  sites  and  also  at  neighboring  sites  in  the  (1,0)  direction. 

\ 

Despite  these  caveats,  we  feel  the  shift  in  the  peak  of  Tj  to  a  lower  0  value 

in  going  from  the  flat  to  stepped  surface  is  compelling  evidence  for  the  (+  +) 

2 

case.  This  conclusion  is  strongly  buttressed  by  the  experimental  results  for 
D=24.  Here  the  peak  is  at  a  higher  0  value  (0  =  0.44),  but  still  less  than  0.5, 
which  is  consistent  with  our  interpretation. 


At  low  coverage 4 j  increases  smoothly  for  either  (+  +)  case.  For  (+  -) 

2 

one  would  expect  it  to  grow  smoothly  («  6  )  for  0  <  1/D  as  one  row  of  adatoms 
forms  at  the  attractive  edge  sites.  For  ( — ),  if  we  assume  the  edge  sites 

fill  up  at  equal  rates  as  0  increases,  there  are  two  possibilities.  For  D=even, 

\ 

the  edges  cancel  in  phase,  so  "j  will  be  quite  small  for  0  <  2/D.  For  D=odd, 

i  2 

one  expects  ^  «  0  . 

Actually,  if  one  or  both  of  the  terrace  sites  were  attractive,  the  low  0 

behavior  could  be  more  complicated^.  Since  the  substrate  structure  is  different 

at  edge  sites  than  elsewhere  on  the  terrace,  the  adatom-adatom  (AA)  interactions 

may  be  different  for  these  sites.  Changes  in  the  AA  energies  for  the  (+  -) 

1  2 

or  ( — )  case  would  not  affect  the  I  domain  LEED  spots  at  low  0  (<  ^  or  £  p, 
respectively).  However  it  is  possible  for  adlayer  LEED  spots  not  observed  on 
the  flat  surface  to  occur,  if  the  AA  energies  impose  a  new  ordering  along  the 
terrace  edges,  or  for  the  II  domain  LEED  spot  to  be  affected.  See  I  for 
further  discussion  of  this  point.  For  this  reason  the  low  6  LEED  behavior 
on  stepped  surfaces  may  be  more  complicated  to  interpret  than  what  occurs  for  0  > 
Calculated  values  for  jj  are  shown  in  Fig.  4.  Note  that  for  the  (+  +) 
cases  we  recover  what  is  observed, =  -  jj  for  low  0  and  .  j  >  j j  at  larger 
coverages.  For  the  ( — )  ((+  -))  case  and  0  <  2/D  (1/0)4^  will  be  very  small 
in  our  model  since  all  the  adatoms  will  be  at  terrace  edge  sites.  (On  the 
real  surface  it  might  be  somewhat  larger  if  there  are  defect  sites  we  have  not 
included  here  that  can  nucleate  a  II  domain).  Thus  the  (+  +)  case  again 
agrees  with  experiment  and  the  others  do  not.  Note  that  none  of  the  cases 
show  the  observed  minimum  in  near  the  maximum  in  ~j.  The  reason  for  this  may  be 
due  to  the  imposition  of  periodic  boundary  conditions  in  the  [T  1  1J  direction 
in  our  model.  Different  boundary  conditions  might  cause  II  domains  to  nucleate 
more  easily  at  low  0,  only  to  disappear  as  0  increases  and  "crowding"  in  a  long 
narrow  terrace  favors  the  I  domain.  However  we  have  not  checked  this  idea  numerically. 


o|ro 
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Now  we  discuss  LEED  spot  splitting.  Although  this  may  be  a  strong  function 

2 

of  8,  experimental  results  were  reported  at  unrecorded  8  values.  We  first 
confine  our  remarks  to  8  =  0.4  since  the  remarks  on  splitting  in  Ref.  2  ap¬ 
parently  refer  to  Fig.  2  of  that  work  which  is  at  this  coverage.  We  found  the 
I  domain  spot  to  be  strongly  split  for  all  (+  -)  and  (-  -)  cases  treated.  If 

,r  1  ~ 

the  spot  center  is  at  k  and  the  maximum  at  k,  =  k_  +  —  yr  aa  measure  of  the 

r  —o  — l  — o  a  D 

intensity  distribution  in  the  split  peaks  is  the  ratio 

R  -  x(k,)  /r  (  k.l  (s) 

For  the  (+  -)  case  we  found  R  =  55  for  D  =  10,  R  =  11  for  D  =  9.  For  the  (+  -) 
case,  R  =  69  for  D  =  10,  R  =  344  for  D  =  9.  For  the  (+  +)  case  R  =  191  for 
D  =  9  while  R  =  0.7  for  D  =  10.  However,  the  real  stepped  surface  must  have  a 
distribution  of  even  and  odd  D,  so  that  the  observed  spot  splitting  is  not  in¬ 
consistent  with  the  (+  +)  value  of  R  for  D  =  10.  For  the  II  domain,  note  that  <p(kj> 
is  the  Fourier  transform  of  <n^>.  Using  periodic  boundary  conditions  this 
average  is  independent  of  the  position  of  site  i  in  the  [T  11J  direction.  Hence 
<p(k)>  =  0  for  the  (1/2,  1/2)  beam  and  for  any  beam  whose  displacement  is  in  the 

(1,1)  reciprocal  space  direction  from  this  point.  So  the  contribution  of 
2  ! 

| <p( k_)> |  to  I(k)  (or  ~jj)  is  very  small  and  one  would  expect  no  splitting  of 

the  II  domain  beam,  as  observed  experimentally.  On  the  real  surface  there  are 

no  periodic  boundary  conditions  but  a  random  distribution  of  terrace  end  locations 
_  2 
in  the  [1  11]  direction  will  have  the  same  effect  on  |<p>|  . 

As  pointed  out  in  I,  the  LEED  spot  splitting  can  be  a  strong  function  of 

8.  This  general  conclusion  is  born  out  by  our  numerical  results.  For  the  I 

domain,  the  ratio  R  in  Eq (17)  for  (+  +),  D  =  10  is  0.95  at  8  =  .05,  0.67  for 

8  =  0.4,  and  1.11  for  8  =  .55.  For  D  =  9,  R  =  0.91  for  8  =  0.03,  191  for  8  =  0.41, 

1577  for  8  =  .45,  and  6  for  8  =  0.63. 


VI.  Prospects  for  Further  Work 

We  have  considered  the  general  theory  of  LEED  scattering  from  certain  types 
of  overlayers  on  stepped  surfaces  including  the  effects  of  disorder  in  the 
overlayer  for  the  first  time.  By  comparison  with  experimental  results  of  0  on  a 
stepped  W(110)  surface,  we  have  demonstrated  the  adsorption  at  terrace  edge  sites  on 
this  particular  surface  is  less  strong  than  elsewhere  on  the  terrace.  This 
conclusion  rests  most  strongly  on  a  single  experimental  fact:  that  the  peak 
intensity  for  LEED  scattering  from  the  I  domain  (see  Section  V  and  Fig.  2) 
occurs  at  a  coverage  0  less  than  that  on  the  flat  surface.  Referring  to  the 
discussion  in  I  and  Section  V,  it  seems  clear  that  the  three  possibilities  for 
terrace  edge  site  adsorption  energy  changes,  (+  +),  (+  -)  and  (-  -)  should  be 
easily  distinguishable  by  a  simple  comparison  of  experimental  results  on  the 
stepped  and  flat  surfaces.  This  indeed  seems  to  be  the  case,  and  we  hope  to 
spell  out  the  full  details  elsewhere.  For  present  purposes  we  confine  our¬ 
selves  to  the  following  remarks,  (i).  The  conclusion  of  Section  V  regarding 
the  behavior  of  -j  near  its  maximum  can  be  obtained  by  a  simple  T  =  0  calculation 
very  similar  to  that  reported  in  I.  This  is  reasonable,  since  the  statistical 
behavior  near  a  complete  overlayer  should  not  be  strongly  influenced  by  entropic 
effects  in  a  finite  system  of  this  kind.  (ii).  The  three  cases  (+  +),  (+  -)  and 
(-  -)  almost,  but  not  quite,  correspond  to  -j  peaking  at  0  =  0.5-jj,  0  =  0.5,  and 
0  =  0.5+^j  respectively.  The  "not  quite"  refers  to  odd-even  effects  (see  the  (+  -) 

cases  in  Fig.  4,  for  instance).  However  this  problem  can  be  overcome  by  considering 
2  2 

< | p |  >  -  | <p> |  ,  a  quantity  which  is  experimentally  accessible,  as  follows  from 
Eq(5).  (iii).  For  lattice  gas  models  with  even  AA  interactions  only,  this  last 
conclusion  also  holds  for  T  >  0  by  an  application  of  up-down  symmetry  in  the 
corresponding  Isinq  model.  If  odd  interactions  are  important,  as  appears  to  be 
generally  the  case^  for  indirect  interactions  in  chemisorption  systems,  up-down 


symmetry  no  longer  obtains.  However  we  expect  the  method  will  still  go  through 


since  the  shift  in  0  dependence  is  really  due  to  stearic  effects,  and  so 
should  have  consequences  only  weakly  dependent  on  interaction  symmetry. 
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Note  added  in  proof 

Recent  work  by  Besocke  and  Wagner^  lends  support  to  our  conclusions 
about  the  change  in  0  binding  energy  on  the  stepped  W(110)  surface  of  ref.  2 
(and  Fig.  2).  Ref.  11  examines  the  coadsorption  of  0^  and  on  two  stepped  W(110) 
surfaces.  One  of  these  has  step  edges  in  the  [1  10]  direction,  the  other  in 
the  [001]  direction.  It  is  concluded  that  N2  adsorbs  dissociate^ly  at  terrace 
edge  sites  in  either  case.  Coadsorption  of  0  poisons  the  N2  dissociation, 
more  strongly  in  the  [T  10]  case  than  the  [001]  case.  This  implies  that  the 
0  binding  energy  for  the  former,  at  edge  sites  on  at  least  one  side  of  the 
terrace,  is  larger  than  on  flat  W(110);  and  that  the  0  binding  energy  is 
probably  slightly  larger  in  the  latter  case.  The  distance  between  W  atoms 
along  a  step  edge  decreases  progressively  as  one  goes  from  the  [T  10]  to  [001] 
to  the  present  (close-pack  direction)  case.  Hence  a  progressive  decrease  in  0 
binding  energy  at  the  terrace  edge  sites  presents  a  consistent  trend. 
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Figure  Captions 

1.  Adatom-Adatom  interaction  energies  for  0/W(110)  as  determined  by  Ching 
et  al^.  Their  values  are  =  -0.072  eV,  Eg  =  +0.08  eV,  E3  =  -0.049  eV. 

2.  The  stepped  W(110)  surface  used  in  Ref.  2.  Large  circles:  W  substrate 
atoms.  Small  circles:  0  adatoms.  Dark  circles  illustrate  overlayer 
domain  I,  cross-hatched  circles  domain  II. 

3.  Experimental  results^  for  integrated  LEED  intensities as  a  function  of 
coverage.  Circles  4^,  crosses^ j. 

4.  Monte  Carlo  results  for  the  integreated  LEED  intensities^/N^  (circles) 

and4^j/NsN2  (crosses)  for  various  terrace  widths  and  terrace  edge  ad¬ 
sorption  energy  changes  (see  text),  (a)  9(++),  (b)  10(++),  (c)  9(— ), 

(d)  1 0( — ) ,  (e)  9(+-) ,  (f)  10(+-) . 
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